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ABSTRACT 

We present new measurements of the cosmic microwave background (CMB) polarization from the final season 
of the Cosmic Anisotropy Polarization MAPper (CAPMAP). The data set was obtained in winter 2004-2005 
with the 7 m antenna in Crawford Hill, New Jersey, from 12 W-band (84-100 GHz) and 4 Q-band (36^5 GHz) 
correlation polarimeters with 3.3' and 6.5' beamsizes, respectively. After selection criteria were applied, 956 
(939) hours of data survived for analysis of W-band (Q-band) data. Two independent and complementary 
pipelines produced results in excellent agreement with each other. A broad suite of null tests as well as extensive 
simulations showed that systematic errors were minimal, and a comparison of the W-band and Q-band sky maps 
revealed no contamination from galactic foregrounds. We report the is-mode and Z?-mode power spectra in 7 
bands in the range 200 < ^ < 3000, extending the range of previous measurements to higher £. The E-mode 
spectrum, which is detected at 1 Ict significance, is in agreement with cosmological predictions and with previous 
work at other frequencies and angular resolutions. The BB power spectrum provides one of the best limits to 
date on fi-mode power at 4.8 /iK^ (95% confidence). 

Subject headings: cosmology: cosmic microwave background — cosmology: observations — polarization 



1. INTRODUCTION 

The CMB has not yet yielded all its secrets, particularly from 
its polarization. Although the outlines of the power spectrum 
of the positive-parity E-modes have emerged since their first 
detection five years ago ( Kovac et al. 2002 1, its details have not, 
and the negative-parity B-modes remain entirely elusive. The 
principle sources of /s-modes are the same perturbations in 
the primeval plasma that gave rise to temperature anisotropics. 
Since the latter are well-measured, robust predictions exist for 
the cosmological EE spectrum on most angular scales. 

To date, EE data from some seven collaborations are in 
agreement with predictions. Nonetheless, more detailed char- 
acterization of the EE spectrum, one of the very few probes 
of the early universe, is important and may reveal surprises. 
Additionally, these efforts lay groundwork for future detection 
of B-modes, which are compelling targets. At the subdegree 
angular scales studied by CAPMAP, B-modes should arise 
from gravitational lensing of /s-modes, and thus provide a 
new channel for examining the intervening energy content of 
the universe. At large angular scales their sole extragalactic 
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sources are gravitational waves from the earliest instants of 
time; detection of those primordial B-modes would give a 
measure of the energy scale of inflation itself. 

Here we report further measurements of the EE spectrum 
and new constraints on the B-mode power from the final season 
of CAPMAP. These data derive from a unique combination of 
observing frequencies and angular resolution, and in particular 
extend to larger multipoles than previous work. CAPMAP 
measured the CMB polarization at 40 GHz and 90 GHz simul- 
taneously, two frequencies chosen to bracket the foreground 
frequency minimum observed by WMAP and to permit new 
constraints on foreground contaminants. Moreover, the mea- 
surements were made at very small angular scales using corre- 
lation polarimeters rather than bolometers or interferometers. 
Thus, we report here data complementary to that from recent 
work (including our own): ,Ade et al. ( 2007|l; Barkats et al 



( 2004 ) ; Leit ch et al.|d2005]l; (Montroy et al.|(|2()()6| ); |Page et al. 



( i2007j ; iSievers et al.| ( |2007l l; |Wu et al.| ( |2007F 

The reported power spectra, drawn from two complementary 
analysis pipelines, are robust, as evidenced by results from 
a suite of 72 null tests on the polarization maps, and from 
extensive simulations of possible systematic effects. 

We begin with description of the experimental configuration 
(^, the CMB observations (Q, and the methods of calibra- 
tion (Q, before moving on to explanation of the data analysis 
(Q. q6]presents the polarization power spectra and the results 
of the null tests, while ^ and ^discuss evaluation of system- 
atic errors and foreground limits. We finish with comments on 
internal consistency (Q followed by conclusions (|T0|. 

2. INSTRUMENT DESCRIPTION 

The CAPMAP instrum ent is described in IBarkats et al.l 
P005| l and summarized in [Barkats et al.| ( r2004[ ). Below we 
outline the final configuration used for the observations re- 
ported herein, highlighting improvements compared to the first 
season. The CAPMAP instrument was optimized to measure 
the EE and BB power spectra, sacrificing sensitivity to TT and 
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Fig. 1. — CAPMAP an'ay configuration and the ring scan. Tiie instantaneous 
locations of tlie 16 radiometers are shown with black (gray) circles for W-band 
(Q-band). The circle diameters show the 3.3' (6.5') beams. The telescope 
moved so that the center of the array continuously traced a constant-declination 
circle about the NCP. Dashed circles indicate sky coverage: each is the path 
followed by one radiometer during a single 21s cycle. Each radiometer 
timestream is considered to be a function of the variable d. 



TE. 

CAPMAP comprised 16 heterodyne correlation polarimeters 
coupled via corrugated feed horns and hig h-density polyet hy- 
lene lenses to the 7 m Cassegrain antenna ( Chu et al.|1978| l in 
Crawford Hill, NJ (W74°ll'll", N40°23'31"). 12 receivers 
operated at W-band (84-100 GHz) with beam FWHM 3.3', 
while the other four operated at Q-band (35^6 GHz) with 
beam FWHM 6.5'. The receivers were organized into four 
groups (numbered to 3), each including one Q-band and three 
W-band receivers (labeled B, C and D) housed in a cryostat, 
which cooled the high-electron-mobility transistor low-noise 
amplifers (LNAs) and feedhorns to 25 K, and the lenses to 
80 K. 

For each radiometer, the circular feed coupled two orthogo- 
nal components of the incident electric field via an orthomode 
transducer into LNAs in each arm. The LNAs, provided by 
JPL, were based on monolithic microwave integrated circuits 
(MMlCs) and had typical noise temperatures of 65 K for W- 
band and 35 K for Q-band. After amplification and bandpass 
filtering, mixers downconverted the radio frequency signals 
to intermediate frequencies (IF) in the range 2-18 GHz with 
local oscillators (LOs) operating at 82 GHz or 30.5 GHz. The 
LO signal for one arm of each polarimeter was phase-switched 
at 4 kHz to suppress 1// noise. The three W-band radiometers 
in each cryostat shared a single LO. 

After IF amplification, the W-band (Q-band) signals were 
divided into three (two) frequency subbands per polarimeter, 
for a total of 44 channels. For each, signals from the two arms 
were combined in an analog multiplier and then preamplified 
before digitization. Thus each polarimeter timestream was 
proportional to the product of the two orthogonal modes of the 
incoming electric field. At a given parallactic angle ?/;, each 
CAPMAP polarimeter measured one linear combination of the 



Stokes parameters Q and U, 

F = 2cos[2(^ + ry)] +f/sin[2('(/' + ?7)], (1) 

where the angle ry, the detector polarization angle, varied for 
each channel but was nominally 45°. 

The subbands were labeled SO, SI and (for W-band) S2. The 
two orthogonally-polarized total power signals were tapped off 
each arm of a polarimeter prior to the frequency subdivision, 
and recorded with detector diodes. Because they were not 
phase-switched, these total power channels were too noisy to 
provide CMB data, but they were useful for assessing radiome- 
ter health, monitoring the atmospheric noise for data selection, 
measuring the properties of the beams, and determining the 
atmospheric opacity. 

The active components (the cold LNAs, the warm LOs, 
the IF amplifiers and the pre-amplifiers) were temperature- 
controlled to stabilize the responses of the radiometers. These 
servo systems did not perform consistently, however In partic- 
ular, their temperature set-points had to be increased multiple 
times, resulting in changes of 10-20 K from the start to the end 
of the observing season. Moreover, the temperatures occasion- 
ally drifted for 5-10 hour stretches in response to large swings 
in the ambient temperature. These drifts were typically at the 
level of 0.3 K hr"' but reached as high as 1-2 K hr"' in some 
cases. As described in ^ these temperature instabilities led 
to a systematic error in the experiment, but it was considerably 
smaller than the statistical error. 

Both the polarization and total power channels were sam- 
pled at 100 kHz by a commercial sigma-delta ADC board, as 
described in Barkats et al. ( 2005) 1. In lieu of recording high- 
speed data as in the first observing season of CAPMAP, the 
polarization channels were digitally demodulated both in phase 
and 90° out of phase with the 4 kHz phase switch clock. The 
out-of-phase (quadrature) data lacked signal from the sky, but 
had the same white noise properties as the in-phase data and 
so were used for systematic checks. A small amount of sig- 
nal (typically < 2%) did leak into the quadrature data, due to 
phase shifts in the data relative to the digitized clock used for 
demodulation, but this led to a negligible sensitivity reduction 
for the in-phase data. The data rate was 2.9 GB per day. 

The array configuration and scan strategy, illustrated in Fig- 
ure [T| were designed to have high symmetry, provide highly 
cross-linked maps, and constrain Q and U tightly with uni- 
form parallactic angle coverage. For CMB observations, the 
telescope's optical axis described a circle around the north ce- 
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2. — Parallactic angle coverage. The solid (dashed) lines show the 
paths of the 12 W-band (4 Q-band) radiometers during a single 21 s cycle. 
Pixels at each declination were observed at a variety of parallactic angles, 
overconstraining the Q and U Stokes parameters everywhere in the map. 
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lestial pole (NCP) with a radius r = 0.725° and a period of 21 s. 
The observed position is parameterized by the ring variable 6 
shown in Figure[T| This scanning method surveyed a 3° disc 
centered on the NCP, covering the entire region every six hours 
with good noise uniformity (±50%) and an even parallactic 
angle distribution (Fig.[2|. The 7 m antenna provided a large 
and uniform focal plane so that even the outermost receivers 
had essentially Gaussian main beam profiles. 

The fact that the elevation varied continuously throughout 
the scan provided several side benefits. The ring motion pro- 
duced a sin 9 modulation in the total power channels with 
amplitude given by ( [McMahon |2006| l 

Ai =-r7;tmCOt(e/o)csc(e/o)rexp[-rcsc(e/o)], (2) 

where r is the opacity of the atmosphere, ratm is its tempera- 
ture, elo = 40.3° is the central elevation of the ring scan, and 
r its radius. The total power channels thus continuously mon- 
itored the opacity and the stability of the atmosphere. Since 
the atmospheric signal is an unpolarized intensity /, sinusoidal 
responses in the polarimetry channels measured the small leak- 
age of / into polarization, which we refer to as / ^ 2 leakage. 
As discussed further in ^5.1| and ^5.2| the atmospheric signal 
can be cleanly removed from the polarization channels (with a 
15% loss of sensitivity). 

We note here the significant improvements in the final in- 
strument configuration for CAPMAP since the first season 
(reported in B arkats et al.| (2004 1). The original four W-band 
receivers have been supplemented by an additional eight W- 
band and four Q-band receivers. Data were collected with the 
ring scan strategy described above, rather than with a constant- 
elevation scan. The lenses were upgraded to reduce I ^ Q 
leakage, primarily by improving their anti-reflection (AR) coat- 
ing. The original coatings comprised birefringent grooves, 
which were replaced with a polarization-symmetric square 
grid of round holes; the resulting improvement is reported in 
^4.3[ Finally, an order-of-magnitude reduction of scan syn- 
chronous signals from telescope sidelobes was accomplished 
in a four-pronged approach (McMahon|2006| l : 1) a co-moving 
ground screen was added between the receivers and the sec- 
ondary mirror; 2) a baffle was added covering the edges of 
the secondary mirror, redirecting rays otherwise headed for 
the ground back onto the primary and thence to the sky; 3) 
gaps between panels of the primary mirror were sealed with 
copper tape, measurably reducing diffractive sidelobes from 
those gaps; and 4) several nearby trees that had previously 
illuminated sidelobes were removed. 

These modifications reduced the linear dependence of the po- 
larimeter response on azimuth by an order of magnitude, from 
500 /^K deg"' to 50 /xK deg"'. Moreover, whereas these slopes 
showed variations of « 200 /iK deg"' day"' during the first 
season, any remaining variations were no longer measurable. 
Concluding remarks on the efficacy of this sidelobe-reduction 
effort are given in ^5.1| and ^ 

3. CMB OBSERVATIONS 

Between 13 December 2004 and 28 April 2005, CAPMAP 
recorded 1658 hours of CMB data: a 53% livetime. Another 
3% of the time was devoted to calibration data. The majority 
of the downtime was due to snow, rain, and heavy clouds, but 
power outages, telescope repair, and other gross failures pre- 
vented observations about 10% of the time. We downselected 
the CMB data as described in ^5.1[ primarily removing more 
periods of bad weather The final sample comprised data from 
29% of the time in the interval above. 



We divided the data into three periods (I, II, and III) of 
roughly equal length and sensitivity, separated by power out- 
ages. These periods are visible in Figures |4] and [5] For the 
purposes of data selection and quality inspection, we further 
subdivided the data into nine subperiods. We also designated 
one long continuous segment of one of the subperiods as the 
reference period (RP), used for initial exploration of the data. 

4. CALIBRATION 

In addition to CMB data, calibration data were taken before, 
during, and after the season. These data measured parameters 
necessary for translating the CMB timestreams into polariza- 
tion power spectra, including the coefficients in the pointing 
model, the beam sizes and shapes, the detector polarization 
angles, and the calibration coefficients C relating the detector 
voltage outputs to thermal CMB fluctuations. These last re- 
quire measurement of both the atmospheric opacity t and the 
detector responsivities TZ (with units V K"') to thermal sources 
below the atmosphere: 

C = 7^/exp[rcsc(e/)], (3) 

where el is the elevation of the observation. We refer to the 
denominator in equation (|3]l as the absorption correction. The 
primary calibration method for TZ required evaluation of tele- 
scope beam efficiencies. We also used the calibration measure- 
ments to determine the level of / — > Q leakage. 

Below we describe the measurements of all these quantities. 
We conclude with a description of the detector responsivity 
measurements, and the model used to extrapolate those mea- 
surements to all times during the observing season. The impact 
of uncertainties in determining all these parameters is assessed 

4.1. Pointing Model 

To determine where the receivers were pointed during CMB 
and calibration observations, we undertook measurements of 
both the global pointing of the array center, and the relative 
offsets of each radiometer from a fixed position in the focal 
plane. 

The global pointing relied on a 10-parameter pointing model 
( Meeks et al.|1968| , constrained by two sets of radio data: one 
obtained on 31 Dec 2003 with 154 source observations, and 
another over 26-27 Nov 2004 with 50 observations. In each 
session, a single Q-band receiver was repeatedly rastered over 
a succession of five sources (Cas A, Cyg A, Tau A, OMC 1 , 
and Jupiter) as the sky rotated. Each source map was fit to a 
Gaussian beam shape, and the recovered center in telescope 
coordinates was compared to the known source position. The 
resulting pointing offsets were used to fit the pointing model 
parameters. The residuals from the best fit had an RMS of 29" 
in elevation and 18" in cross-elevation. The results from both 
data sets were statistically consistent, showing the pointing 
solution to be stable in time. 

The relative offsets of the radiometers were determined from 
1 1 observations on 26-27 Nov 2004, in which the full array 
was rastered across Jupiter. Again, source maps were fit to 
Gaussian beams. The measured offsets agreed well with op- 
tical simulations, to within the known alignment tolerances. 
The uncertainties in the offsets were typically 4" and always 
smaller than 7". 

4.2. Beam Measurements 

The beam size and shape for each radiometer were measured 
with observations of Jupiter in the total power channels. Two 
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scan strategies were used, both optimized for efficiency by 
concentrating integration time on single receivers. The first 
strategy consisted of a tight raster scan with a 21' width and 
a 0.3' elevation step. The second, a point and integrate scan, 
used discrete steps spaced on a 1.5' grid with a 2 s integration 
time at each point, bracketed in time by 2 s off-source mea- 
surements. For both scan types, the resulting maps were fit to 
two-dimensional Gaussians to extract the beam parameters. 

The two scan strategies produced consistent estimates of the 
mean beam FWHM in each frequency channel: 3.3' and 6.5' 
for W- and Q-band respectively. The standard deviation of the 
beam size measurements of all W-band (Q-band) radiometers 
was 2% (3%); all were within 5% of the mean in each case. We 
note that the 2Q beam — the highest in the focal plane — was 
slightly truncated as it left the vertex cab, which elongated it 
and increased its beam size by about 6% compared to those 
from the other three Q radiometers. 

Since the center frequencies of the polarimeter subbands 
differed from those of the broadband total power channels, and 
since the optics were near diffraction-limited, the beam sizes 
of the high and low subbands differed by about ±2% from the 
measured total power values ( |McMahon|2006| l . 

The beams were all quite close to circular, with elongations 
less than 8%. Elongation is defined as (b-a)/{b+a), where a 
and b are the beam radii along the principal axes. 

4.3. I Q Leakage 

Jupiter measurements also constrained the I ^ Q leakage 
parameters of each receiver; the point and integrate scan was 
optimized specifically for this purpose. Since the 1 // noise of 
the polarization channels was negligible over the four seconds 
per point in the map, it was possible to integrate down to 
measure the very small polarized signals due to optical and 
radiometric I ^ Q. These measurements showed that the 
m onopole, dipole, and q uadrupole terms (see |Hu et al.| ( |2003 i 
or |Barkats et al.| ( [2005) 1 for the definitions) were below -23 
dB, -15 dB, and -20 dB respectively. These levels are small 
enough to neglect, as discussed in ^ The new lens AR 
coatings described in ^improved the quadrupole suppression 
by a factor of ten compared to the original grooved coatings 
used during CAPMAP's first season. 

The monopole leakage terms were also measured by com- 
paring each polarimeter' s response to the atmospheric sinf? 
variation in the CMB ring scan to the response from the total 
power channels, as mentioned in ^ These measurements 
were consistent with those from Jupiter. 

4.4. Detector Polarization Angle Measurements 

Observations of the moon were used to determine the detec- 
tor polarization angles 77. Thermal radiation from the moon is 
polarized according to its p ropa gation direction relative to the 
normal to the surface (see |4.7| for a related discussion). Thus, 
the polarization pattern of the moon is oriented radially and 
peaks in amplitude near the edge of the moon's disc. Projected 
onto the sensitivity axes of a CAPMAP polarimeter, the radial 
pattern was seen as a quadrupole, which we fitted to find the 
angle 77 as illustrated in Figure [3] The measurement error on 77 
ranged from 1° to 2.5°; the resulting systematic uncertainty in 
the final power spectra was small, as discussed in ^ 

These measurements of 77 also revealed that one radiometer 
was installed nearly 20° rotated from its design orientation, 
which was confirmed during decommissioning. This affected 
the parallactic angle coverage but was taken into account in 
the analysis. 
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Fig. 3. — Polarized map of the moon from chamiel 2BS0. The dashed 
contours indicate negative signal levels. The polarization angle for this channel 
was fitted to be 44.7° it 1.5°. The moon was not full, as can be seen by 
the envelope on the quadrupole pattern, which features deeper lobes on the 
left side; this effect was taken into account when determining the detector 
polarization angles. 

Observations of Tau A at multiple parallactic angles were 
also fitted for 77 and gave consistent results. 

4.5. Opacity Determination 

The ring scan permitted continuous monitoring of the opac- 
ity during CMB observations. Ten times during the season, 
the opacity was independently measured from Jupiter scans 
followed by sky dips, in which the telescope was scanned in 
elevation from the zenith to near the horizon and back. The 
combination of the Jupiter and sky dip data were analyzed to 
yield both total power responsivity estimates and opacity esti- 
mates. These opacity measurements confirmed the continuous 
monitoring results, as shown in Figure [4] 

The clear sky brightness temperature data from the Geosta- 
tionary Operational Environmental Satellite'^ were used to 
estimate Tatm, which varied from 255-295 K over the course 
of the season. The amplitude A 1 (eq. (|2]i) was averaged each 
hour across the 7 (6) most stable W-band (Q-band) total power 
channels. Numerical evaluation of equation (|2]) resulted in the 
hourly estimates of r plotted in Figure [4j which were splined 
to the 21s intervals on which the calibration was applied. The 
average value for r after data selection was 0.135 (0.092) for 
W-band (Q-band), corresponding to an absorption correction 
of 1.23 (1.15). 

Measurements of r were limited to 10% accuracy by uncer- 
tainties in the total power responsivities used to convert A 1 to 
temperature units. The resulting uncertainty in the absorption 
correction depended on r, but was in the worst case 5% (2%) 
for W-band (Q-band). The season average was a 1.7% (1.6%) 
uncertainty in the absolute calibration for W-band (Q-band); 
the impact of this uncertainty is discussed in ^ 

These t measurements agreed with the Jupiter plus sky dip 
measurements at the 5% (7%) RMS level for W-band (Q-band), 
consistent with errors in the latter 

4.6. Detector Responsivity Measurements 

As in the CAPMAP first season, the primary measurement 
of the responsivity 7?, in V K"' for each polarimeter was made 
using a nutating aluminum flat (the chopper plate) installed in 
front of the secondary mirror ( Barkats et al. 2004 2005) . When 
radiation from the sky is reflected obUquely by the plate, the 
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Fig. 4. — The atmospheric opacity (left axis) and the con'esponding absorption correction (right axis) for all retained data as a function of time since the start 
of the observing season. Shown are the separate W- and Q-band measurements, from both the radiometer data (labeled sin 9) and the 10 Jupiter plus sky dip 
observations. Shadings indicate the three periods and the reference period discussed in ^ 



non-zero resistivity of aluminum leads to different reflection 
coefficients for the transverse electric and transverse magnetic 
boundary conditions, producing a polarized signal from ini- 
tially unpolarized radiation. The same effect also causes the 
thermal radiation emitted from the plate to be polarized; the 
net signal incident on a polarimeter is given in temperature 
units by 



P(a): 



4^7reo/9J^F(7;ky - Tpi^ts)a, (4) 



where a is the nutation angle of the plate, v is the central 
(band-averaged) frequency, p is the resistivity of the plate, Y 
is a geometrical factor of order unity, Tpiate is the temperature 
of the aluminum, and T^ky is the sum of the atmospheric tem- 
perature and the CMB temperature, Tcmb- The quantity in 
brackets involving x = hv/kTcuE is the thermodynamic cor- 
rection, which puts the responsivity into temperature units 
appropriate for fluctuations in a blackbody source at Tcmb- 
The full nutation range is Aa = 13°, producing a modulated 
signal with peak-to-peak amplitude AP approximately 90 mK 
for typical sky and plate temperatures. The responsivity is 
found from TZ = e,„Ay/Af , where e,„ is the beam efficiency 
( Kraus|1966) l, described in the following paragraphs, and AV 
is the measured amplitude from the polarimeter in volts. The 
uncertainty on AP is 7.5%, dominated by a 5% uncertainty 
on the value of y^, and 4% uncertainties on each of Aa and 

(7sky~ Tplate)- 

Because the chopper plate intercepted the beam before the 
primary and secondary mirrors, each responsivity measure- 
ment from it had to be corrected to account for the telescope 
beam efficiency e„j. A random surface error RMS of as de- 
creases the responsivity to a beam-filling thermal source by 
the Ruze factor ( Ruze|1966) l at the wavelength A: 



: exp 



(5) 



We estimated e„, and an associated uncertainty for each fre- 
quency channel using simulations constrained by measure- 
ments. The measurements included observations of the moon 
with the CAPMAP receivers, and original surface measure- 
ments from jiaret^([T978]l. Specifically, |Chu et al.| ( |l978| ) 
reported the range of measured RMS on the individual panels 
of the primary mirror before installation (ap) and the existence 
of a 150-200 jiva step {d) between the top and bottom halves 
of the mirror. Furthermore, they measured e,„ at 100 GHz with 
a prime-focus feed system and inferred cr, = 100 fiva with a 
3(7 error of 30 fjm. The secondary mirror surface errors were 
negligible. 



The simulations treated the phase error on a 2048 x 2048 
grid covering the telescope aperture. Each simulation drew ap 
for each panel from a Gaussian distribution = 40 ± 20 /im, 
and d from 175 ±25 /im. Moreover, each panel was consid- 
ered to have an additional fourth-order polynomial distortion, 
consistent with its having an adjuster in each comer; these dis- 
tortions were constrained by the requirement that the overall 
surface RMS for all the gridded points be drawn from a flat 
distribution between 70-130 fim. In ^we describe a method, 
using the sidelobes from these simulations, to estimate the 
systematic effect from ignoring them in the analysis. 

The simulations produced beam patterns consistent with the 
100-GHz one plotted in Chu et al. ( 1978| l when the appropriate 
feed pattern was used; that feed resulted in a 3 dB illumination 
contour that encompassed the central eight of the 27 panels 
on the mirror. By contrast, the CAPMAP W-band receivers 
each illuminated approximately two panels at that level. For 
each W-band receiver, many realizations of the simulated beam 
patterns were each convolved with a 0.5° thermal source and 
compared to the moon data, which comprised measurements of 
the total power response out to 1° from the moon center. Good 
qualitative agreement was obtained, but only for simulations 
which conformed to the requirement that it, was contained in 
[70, 130] /im. From consideration of simulations of all the 
W-band receivers, we estimated e,„ = 0.89 ±0.05. We note that 
this estimate is in good agreement with the value 0.87 ± 0.03 
obtained from inserting as= 100 ± 10 /im into equation (|5|. 

For Q-band, similar simulations were made but studied in 
less detail since the final estimated beam efficiency was near 
unity: e,„ = 0.97 ±0.01. The final absolute calibration errors 
from combining the AP and e,„ uncertainties in quadrature 
were 9.4% (7.6%) for W-band (Q-band). 

The responsivity for each of the 44 polarization channels 
was determined from chopper plate observations at 10 separate 
times, spaced approximately evenly throughout the observing 
season. In each of the 10 measurements, the chopper plate was 
run for approximately 30 min and the amplitude was fit in 6 s 
intervals containing w 3 nutation periods each, determining 
the responsivities with better than 1 % statistical precision. 

The responsivity of each channel was also measured through 
6-7 observations of the polarized source Tau A during and after 
the season. The absolute accuracy of these measurements was 
limited to 25% (14%) in W-band (Q-band) by the uncertainty 
in the polarized flux from Tau A, measured most accurately 
by WMAP ( [Page et al.|2007[ ). The Tau A measurements were 
noisier than the chopper plate results and were made at dif- 
ferent times, but they did confirm the long timescale features 
discussed in 34.71 
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Fig. 5. — Responsivity variations during the observing season. Measurements for two polarization channels from the 10 chopper plate observations are plotted 
as plus symbols, with the < 1% statistical error bars omitted. As discussed in the text, 30% of the channels exhibited a significant responsivit y incr ease during 
period II, which is clearly visible for the channel on the right. The straight lines show the three-parameter baseline model (model 2 described in KtJ, which was 
used to interpolate the responsivity over the full season. The circles show the alternative model 1, in which the variations (other than the period II discrepancy) 
were predicted from measured temperature dependences. The model 1 points have been averaged into large bins during contiguous CMB observation runs, with the 
error bai's indicating the spread in the predicted responsivity due to temperature variations during each run. 

which the worst case was 80%, were related to the power out- 
ages at the boundaries of period 11. The sensitivity of each af- 
fected radiometer in mKy^ did not vary measurably across the 
periods — the shift was only a change in the V K"' responsivity. 
Finally, one radiometer, ID, showed a season-long responsiv- 
ity decrease in all three of its subbands three times larger than 
predicted from its temperature coefficients measured at the end 
of the season. The radiometer otherwise performed normally; 



As described in ^4.2| the beam from radiometer 2Q was 
slightly truncated by the receiver cabin window, and so did not 
illuminate the chopper plate uniformly, which led to overesti- 
mation of the 2Q responsivities (since Y in equation Q was 
decreased by an unknown amount). In fact, the 2Q responsivi- 
ties inferred from the chopper plate measurements implied a 
much better sensitivity than measurements of the LNA noise 
temperatures allowed. Thus, the responsivities for 2QS0 and 
2QS1 were adjusted to match the Tau A measurements. 

Since the temperatures of the active polarimeter components 
were known to have varied throughout the observing season, 
additional chopper plate observations were undertaken at the 
end of the season to measure the linear variation of the re- 
sponsivity of each polarimetry channel to each of the three 
germane servoed temperatures described in ^ In each case, 
the responsivity decreased with increasing temperature. The 
temperature coefficients varied for different channels but were 
typically 1 %/K for the LNA and IF ampUfier temperatures 
and 0.5 %/K for the LO temperatures. 

4.7. Responsivity Model 

The chopper plate measurements revealed that the detec- 
tor responsivities were not constant with time, with typical 
channels varying by 20% relative to their initial values, and 
the worst channels varying by as much as 60%. We describe 
here the time-dependent detector responsivity model used for 
calibration of the CAPMAP timestreams. 

Figure [5] shows the 10 chopper plate measurements for two 
representative channels. Since the temperatures of the active 
polarimeter components were seen to increase throughout the 
season, the approximately steady decrease of the responsivities 
with time was not unexpected. In fact, considering the effects 
of all the temperature coefficients led to a typical prediction of 
a 20% change across the full season, in good agreement with 
observations for most (w 70%) channels. However, complete 
recovery of the responsivity variations on shorter timescales 
was prevented both by occasional periods of grounding prob- 
lems in the readout thermometry, and by phase lags due to the 
poor placement of some thermometers relative to the active 
components. 

The panel on the right in Figure |5]demonstrates an anoma- 
lous shift in the measured responsivity during period II. Such 
shifts affected about a quarter of the channels at the 10% level 
or greater and were not forecast. We suspect the shifts, for 



this decrease was also only a change in the V K" responsivity. 

To properly track the time dependence of the responsivities 
of the channels throughout the run (and to gauge our uncer- 
tainties in the understanding of said time dependence), we 
evaluated two different models for fitting the 10 chopper plate 
measurements and interpolating to each 21s ring cycle. 

For each channel, model I's fit used the three measured 
temperature coefficients (B,) and known housekeeping tem- 
peratures (I/(f)) to constrain two parameters — A, an overall 
constant, and A//, a possible offset for period II only: 



Ri(t)=A+A„ + Y,Bim. 



(6) 



1=1 



Model 2, now with three parameters per channel, replaced 
the last term in equation (|6]l with a simple linear term: 



R2{t) = C+Cu + Dt, 



(7) 



where D is a single negative slope (again different for each 
channel), characterizing the overall responsivity decrease with 
time, and C and C// are constants analogous to A and A// in 
model 1. 

Both models were similar in that they accounted for the pe- 
riod II anomaly and for the gradual decrease of the responsivity 
throughout the observing season. Because model 1 applied the 
temperature coefficients to each ring cycle, it also addressed 
shorter-timescale changes in the responsivities. However, as 
mentioned above, use of the temperature coefficients did not 
faithfully reproduce those changes in all instances. For this 
reason, we used the simplified model 2, illustrated for two 
channels in Figure |5] as our baseline model. The model ac- 
counted for the main features of the responsivity variations, 
with an RMS of 6.5% for the 440 residuals from 44 channels 
and 10 chopper plate measurements. The systematic implica- 
tions of its uncertainties are discussed in ^71 
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5. DATA ANALYSIS 

This section describes the steps in the data analysis culminat- 
ing in polarization power spectra. We discuss first how the data 
were selected for processing, then their reduction into maps, 
and the subsequent estimation of power spectra. Two separate 
analysis pipelines (arbitrarily named 1 and 2) were developed; 
at each step we discuss the approach of each. Throughout 
the analysis, we relied on simulations and consistency checks 
to evaluate data quality, immunity to systematic effects, and 
internal consistency of the pipelines. We conclude the sec- 
tion with descriptions of the simulation tools and the null test 
methodology. 

5.1. Data Selection 

Since CAPMAP collected data from near sea level in the 
Garden State, it was critical to cull data contaminated by large 
spatial and temporal fluctuations in the atmosphere during poor 
weather. In this subsection, we describe all cuts made on the 
data and conclude with some observations about properties of 
the final data set. 

We list here several preliminary cuts, indicating parenthet- 
ically the fraction of the 1658 hours of CMB data discarded 
by each. Data selection began by dividing the raw timestreams 
from each channel into 21s ring scan cycles. Only complete 
cycles were selected for inclusion in the final data set. Cycles 
exhibiting irregular telescope motion were discarded (3.8%), 
as were those during periods when the cryostats were warm 
or thermally unstable (6.9%), or when known electrical or 
receiver malfunctions existed (2.4%). Only for the Q-band 
receivers, the sun passing through narrow, constant far side- 
lobes at specific times appeared as a several-mK signal in the 
polarimeters; cycles were cut according to the solar position 
(7. 1% of the Q-band data only). 

The remaining selection criteria were all derived from the 
receiver data, and were initially developed independently for 
the two pipelines. Each used different estimators of the data 
quality. In cutting, pipeline 1 rejected all polarimeters from 
individual cycles, while pipeline 2 made cuts on contiguous 
sets of 1 1 cycles and allowed for channel-specific cuts. In 
the end, both pipelines converged to a common technique 
described below: pipeline 1 's quality estimator applied on both 
timescales, with channel-specific cuts where appropriate. In 
distinguishing among differing techniques, we were guided by 
the results from the null test suite described in ^5.5| and ^ not 
by the power spectra of the CMB maps. 

To further characterize each cycle, for each total power (TP) 
or polarimetry channel we averaged the data into Ng angular 
bins and fitted to the following five-parameter model: 

M(0) =Ao+Ai sin 9 +A2 cos 6 + A3 sin26' +A4 cos 29. (8) 

For data selection, we used A^g = 18. (In subsequent analysis 
steps, the data were binned into 200 or more angular bins, 
depending on the pipeline.) 

The first two terms in this Fourier mode expansion dom- 
inated the response of every channel: Aq describes the TP 
system temperature (or the polarimeter offset), and Ai the 
channel's first-order response to elevation changes, and thus 
to Tatm as in equation (|2]). Their values for TP channels were 
typically 100 K and 500 mK respectively; for the polarization 
channels they were < 150 mK and < 5 mK, governed by the 
size of the / — > 2 leakage. For a single cycle, the remaining 
three terms were typically consistent with zero. Of these, A2 
and A4 characterize an azimuthal dependence, while A3 is the 
second-order term in the elevation dependence. 
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Fig. 6. — Histogram of the cut variable TP-x^ defined in tlie text. Tlie 
black histogram shows the entire RP (24,000 contiguous 21s cycles). A 
10,000-cycle best-weather subset (blue histogram) was identified by eye from 
the timestream shown in Fig. [7] The distribution drawn in red was obtained 
from the entire RP distribution after applying just the 11-cycle cut (see text); 
it approximates the best-weather distiibution. The single-cycle cut (vertical 
dashed line) removed the long tail. 
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Fig. 7. — Timestream of the TP-X" cut variable defined in the text, during 
the RP. Cycles eliminated by the two-timescale weather cut described in the 
text are indicated in red. The RP was chosen to be a particularly quiet period 
with no gaps, so a smaller fraction of cycles was cut than was typical. The 
dashed lines mark the best-weather subset. 



All subsequent cuts were based on these fits for each cycle, 
most often on their goodness-of-fit statistics. The primary 
cut variable used to reject bad weather was obtained from 
values from a subset of the TP channels, comprising the nine 
W-band channels with the lowest 1// noise, constrained so 
that at least one such channel was in every cryostat. (Since 
all TP channels had 1 // noise in excess of white noise, their 
mean x^ values were all in excess of 13, the number of degrees 
of freedom.) For each subperiod, the nine x^ distributions 
were renormalized to have their peaks at the same values, and 
then the renormalized x^ statistics were averaged to provide 
the cut variable for each cycle, a variable we call the TP-x^- 
Atmospheric contamination in the Q- and W-band TP channels 
was sufficiently correlated that the TP-x^ also worked well for 
Q-band selection. 

The RP (see Q was a particularly clean six-day period of 
data identified early in the analysis and studied extensively for 
preliminary work developing the TP-x^ as a cutting variable. 
For reference, the RP is indicated in Figure |4] Figure |6]plots 
the distribution of the TP-x^ during the RP, while Figure [t] 
plots its timestream and gives an indication of its power for 
data selection. 

As indicated previously, selection on two timescales sharp- 
ened the effectiveness of the cuts. The TP-x^ distribution 
in conjunction with results from the null test suite (^ was 
studied with many combinations of parameters defining the 
single-cycle cut and the longer timescale cuts. The chosen 
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Fig. 8. — The distribution of the five-parameter x" statistic is shown in the upper panels for cycles from one of the nine subperiods, with about 80 hours of data, 
for two polarimeter channels after selection cuts. A distribution with 13 degrees of freedom is overlaid as a smooth curve. The channel on the left, a typical one, 
has a mean x~ of 12.95, consistent with the expected value 12.94 it 0.04 (for a distribution truncated at 3(t), while the one on the right, one of the worst channels, 
has a mean of 13.62, which is discrepant at more than 16(j. The lower panels show the noise power spectra for the same two channels. The solid line is a 
three-parameter fit: the white noise level, the knee frequency, and the slope of the 1// component. The dashed vertical lines indicate the knee frequency of each 
channel (4.3 and 13.6 mHz for the left and right panels, respectively); the ring scan frequency can be clearly seen as a high point at roughly 50 mHz. Higher levels 
of 1// noise in the channel on the right increase the variance at the ring scan frequency by 12% over the variance of 100 Hz samples, inflating the x^ statistic. Of 
44 data channels, only 10 had knee frequencies greater than 10 mHz. Although the five-parameter x^ statistic is extraordinarily sensitive to 1// noise, systematic 
studies showed the effect of 1 // noise on the full-season results to be negligible (^. 



parameter set was: any 11 cycles were removed if at least 
seven of them had a TP-x^ > 1-5; and any single cycle was cut 
if its TP-x^ exceeded 2. These cuts were applied uniformly to 
all polarimetry channels, discarding 19.4% of the full data set. 

The five-ring-mode fits to the polarization channels were 
also used for the final cuts on the data.''' For any candidate 
cycle, the sum of the values for all functioning W-band (and, 
separately, Q-band) polarimeters was formed, and if either of 
these sums was > 3a high, that cycle was cut for all polarime- 
ters. Similar cuts were made on the maximum residual to the 
five-parameter fit for all functioning W-band (and separately 
Q-band) polarimeters in a given cycle. For these two cuts, the 
idea was to use the overall goodness of fit for the ensemble of 
polarimeters in a given cycle to keep either all or none of them. 
A further 1.6% of the data was so removed. The next cut was 
based upon the fact that the noise of a properly functioning 
polarimetry channel should depend linearly on the total power 
level. This relationship was examined for each polarimeter, 
for each of the nine subperiods. The linear relationships were 
dominant, but for some isolated periods, cycles were found to 
deviate significantly from the expected behavior As these de- 
viations indicated receiver malfunctioning, those cycles were 
discarded (3.5%). A final cut was invoked to remove glitches 
by eliminating cycles for individual polarimeters with a resid- 
ual deviant at more than 3(7, removing 1 .7% of the Q-band and 
0.3% of the W-band data. 

We close this subsection with some observations about the 
data quality of the polarimeters derived from examination of 
their five-ring-mode fits. In particular, we define two measures 
that were not used in data selection but were used to classify 
data in the null test suite (^. 

In contrast to the TP channels, the distributions from 

We used simulations to verify that the presence of the CMB signal in the 
timestreams did not in any way bias the selection of cycles. 



the fits to most polarimeter channels did closely follow the 
expected \^ distribution for 13 degrees of freedom, verifying 
at the same time our noise model (white noise) and that five 
parameters described the data well. One such channel's distri- 
bution is shown in the top panel on the left in Figure |8] A few 
channels, however — including the one depicted on the right 
in Figure [s] — exhibited distributions shifted slightly high. 
Given the large number of cycles entering these distributions, 
even a shift of less than 1.0 was significant at more than 
10 standard deviations. It was found, not surprisingly, that 
such channels were the ones that showed the worst (highest) 
1 // knees; and that such channels were the same in every one 
of the 9 subperiods. For comparison, Figure|8]shows a noise 
spectrum for each of the two channels in the lower panels. We 
defined a polarimeter' s grand as the mean value of its 
distribution; we tested for contamination associated with high 
grand in ^ We also simulated the effects of 1 // noise 
directly, as described in ^ 

For each polarimeter channel, and again for each subperiod, 
we constructed the 18-element ground-synchronous structure 
vector from an appropriate average of the residuals in each 
6 bin of the ring scan. We evaluated the for the model of 
zero ground-synchronous structure for each polarimeter; the 
W-band channels were consistent with no such signal, but two 
of the eight Q-band channels were not. Those two channels 
showed no evidence, however, for time-variability in their 
ground-synchronous structure. Scan-synchronous pickup from 
telescope sidelobes would be manifest as nonzero ground- 
synchronous structure, so these results are further evidence 
that the telescope sidelobe reduction efforts (^ were largely 
successful. Nonetheless, as described in 5 5.2 the ground- 



synchronous structure mode was projected out of the data 
before power spectrum analysis. 
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5.2. Mapmaking 

The timestream of a given CAPMAP channel can be decom- 
posed into signal plus noise as follows: 

d = Ps + n. (9) 

Here, d denotes the timestream values, represented as a vector 
of length A^tod- The celestial signal in the map space is repre- 
sented as a vector s with 2A^pix components Sp, where the index 
p denotes a combination of a pixel and a choice of Stokes 
parameter The A^(od x 2A^pix pointing matrix is denoted P. The 
coordinate system in which s is defined and the corresponding 
form of the pointing matrix were different for the two pipelines, 
as discussed below. 

The objective of mapmaking is to invert equation (|9| statis- 
tically, arriving at an estimator x for the CMB signal s, with 
associated noise covariance matrix A^. The estimated signal x 
is called a map vector. It is well known that there is an optimal 
choice of unbiased estimator (e.g., Tegmark ( 1997 1), defined 
by the pair of equations 



A^-'x = P^C"'d 



(10) 



assuming that the detector noise n can be modeled as Gaussian 
with covariance matrix C = (nn^). 

The covariance matrix C needed in equation ( [TO] i was de- 
rived from a noise model that assumed white noise. Then, 
certain modes in the data that were to be marginalized in the 
mapmaking were assigned large variance as follows: 



C = Co + e-^R^R + f-'G^G, 



(11) 



where Co was diagonal, e was a small regulating constant, and 
each row of the matrices R, G contained one timestream mode 
to be projected out. The assumption that Co was diagonal — 
that is, that the receiver noise was uncorrelated at different 
times — was verified by data quality studies. The white noise 
level used in Co was independently estimated in each cycle to 
allow for time dependence. The noise estimation method was 
tested extensively in simulations and found to reconstruct the 
true white noise level with acceptable accuracy. The matrix 
R contained ring modes: the lowest five Fourier modes in the 
ring variable 9 in each cycle. The matrix G contained ground- 
synchronous structure modes: one mode for each of Ne bins 
around the ring scan. 

Assigning a large variance to the ring and ground- 
synchronous modes was merely a device for projecting them 
out when making A^~' -weighted maps (eq. ([TOjl). Those modes 
would be removed completely in the limit e — 0. It was prefer- 
able to keep the noise matrix invertible, however, by assigning 
a small nonzero value to e, chosen to be large enough to avoid 
numerical difficulties but still insignificant compared to typical 
elements of C. We verified that changing this parameter over 
three orders of magnitude did not affect the results of either 
pipeline. 

Ring mode removal was necessary to project out the con- 
taminating atmosphere signal from I ^ Q leakage, but it also 
served to highpass filter the timestream, so that the remaining 
noise could be treated as white. The mode removal led to an in- 
crease of the statistical error on the final result in the full range 
2<l< 3000 by 15%. Most I bands were affected only by 
5-10%, but the sensitivity of the lowest band, 2<i< 500, was 
reduced by 70%, and the sensitivity for i < 200 was essentially 
eliminated. 



The noise was verified to be uncorrelated between different 
channels, including different subbands of each polarimeter, 
and consequently we estimated a net data vector x for a given 
subset of W- or Q-band channels by summing the individual 
noise-weighted maps (eq. pO] )) and adding the inverse noise 
covariance matrices. We also combined W- and Q-band chan- 
nels in an enlarged pixel space (2+VK analysis) as described in 

In pipeline 1, the pixel space was a standard pixelization 
in declination and right ascension (RA), using HEALPix'"^ 
with A^side = 2048. At this resolution there were w 10^ pix- 
els in our field, each of diameter 1 .7' or roughly one-half the 
W-band beam FWHM. Mapmaking by brute force applica- 
tion of equation ( [TO] l would require 0{N^^^) matrix operations 
and therefore be prohibitively expensive. However, we found 
a mathematically equivalent algorithm (Smithn2007j ), based 
on repeated application of the Sherman- Woodbury formula 
( Bartlett|1951| l, that was computationally efficient. 

In pipeline 2, the underlying pixel spac e was an azimuthal 
pixelization about the NCP as described by 'Crittenden & Turok 



( 1998 1, with Ns = 90 bins in declination (S) and A^^ = 5 12 bins 
in RA. Then, the largest pixels were at the outer edge of the 
survey region, with diameter 1 .2' or one-third of the W-band 
beam. However, estimation of x was obtained in the half- 
Fourier space ((5, A;ra), where here and below the symbol A:,, 
represents the Fourier transform variable for a coordinate v. 
Prior to construction of the pointing matrix P, the timestream 
data from each polarimetry channel were binned into (9, LST) 
coordinates. Data at fixed 9 were Fourier transformed in LST, 
resulting in a receiver map in (6',feLST)i used as d in equa- 
tion ( [TO] l. Each 9 corresponded to a fixed declination, hour 
angle, and parallactic angle, so ^ra only differed from ^lst by 
a fixed phase for each 9, and P was straightforward to compute. 
Ring mode removal was performed as in pipeline 1, using equa- 
tion ( [TT| , though without the term G^G. Ground-synchronous 
structure corresponds to an LST-independent signal, which 
was removed in these coordinates by discarding the DC mode 

^RA = 0. 

The use of the Fourier transform greatly accelerated the map- 
making process, with evaluation scaling as 0{N^Nj) rather 
than 0{N^Ng). This acceleration came about because in this 
basis the noise covariance matrix was block diagonal — it did 
not connect pixels with different values of A:ra — in the limit 
that the noise in each channel and each 9 bin was independent 
of LST. In fact the noise was only LST-invariant at the w 10% 
level, but extensive simulations showed that the approximation 
of invariant noise resulted in a negligible systematic error, as 
discussed in ^ The half -Fourier coordinate system also diag- 
onalized the theory covariance matrix, yielding a similar speed 
increase for the power spectrum analysis. The power spectrum 
estimation thus broke down into separate subproblems, 
one for each unique value of A;ra. Further acceleration by a 
factor of three was obtained by discarding Fourier modes with 
^RA > 80, which only probe angular scales beyond i = 3000. 



5.3. Power Spectrum Analysis 

The power spectrum analysis was performed by calculating 
the likelihood C, 



-2\n£ = \ndet(N+S) + xHN+Sy^x 



(12) 
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with respect to the theory encoded by 5'. The theory covariance 
matrix 5 can be written 



(13) 



where a, represent M parameters to be estimated from the 
data and the theory basis matrices ^, are calculated from trial 
EE and BB (and sometimes EB) power spectra using standard 
formulas (Zaldarriaga 1998). Results for the parameters a,- are 
reported in ^for several different choices of the number of 
bands M and the trial power spectra. 

To estimate the overall sensitivity and to evaluate the agree- 
m ent of the data with t he concordance model for EE defined 
in Spergel et al. '(2007) and evaluated using CAMB, we per- 
formed two-band fits with one EE and one BB band in the 
range ^ G [2, 3000]. For these, the EE parameter was a dimen- 
sionless multiplier to the concordance model, while the BB 
parameter was a flat bandpower in /iK^. For evaluating null 
maps (see 55.5 and we made similar two-band fits, but 



with both the EE and BB parameters taken as flat bandpowers. 

To estimate spectral features for the final analysis, we formed 
seven bins in £ with width A£ = 300, except we expanded 
the lowest bin {2 < £ < 500) to account for mode removal, 
and we extended the highest bin out to ^ = 3000 to provide 
a check for unexpected power (perhaps from point sources) 
on small angular scales. The choice A£ = 300 resulted in 
small correlations between adjacent bandpowers (typically 5%, 
and 20% for the worst case); the correlations rose sharply for 
narrower bins because of mode-coupling from the finite size 
of the survey region. We estimated EE and BB for each bin as 
flat bandpowers, resulting in 14 bands in total. 

For the Q+W analysis, the pixel space was expanded by 
a factor of two, so that each map pixel appeared once for 
each of the two different beam sizes. The noise was taken, 
as usual, to be uncorrelated among all polarimetry channels, 
but the theory matrix was constrained to represent the same 
underlying CMB sky at both frequencies (properly accounting 
for the beam window function for cross-correlating two beams 
with Gaussian widths (Ji,<T2'- Biiui , ca) = exp 
(tI)/2\). 

From the M-dimensional joint likelihood function 
£(q!i , . . . jttM), we calculated the M one-dimensional 
marginalized functions £,(q!,), from which all statistical 
inferences followed. Each result is quoted as the maximum 
likelihood value with an asymmetric la confidence interval 
calculated as the 68% interval of highest density. For BB 
bands we also calculated 95% confidence upper limits; the 
limits were derived using only the portion of the likelihood 
curves in the physical region of positive bandpower 

The calculation of the joint likelihood was accelerated by 
the method of signal-to-noise eigenmodes ( Bond et al.|1998 1, 
which transforms equation ( [T2] i into a simple algebraic ex- 
pression that determines the likelihood on any dense one- 
dimensional subset of the M-dimensional parameter space with 
only two expensive matrix operations. For the two-band anal- 
ysis, we calculated the joint likelihood on a 60-element grid 
of such subspaces and obtained the marginalized likelihoods 
by numerical integration. For larger M, a useful procedure 
was to use the signal-to-noise eigenmode method to obtain 
M conditional likelihoods £(Q;, |ajy,), in which the likelihood 
for the ith band was calculated as a one-dimensional function 



with the other M- \ bands fixed to the maximum hkelihood 
location. These were only approximations to the true marginal 
distributions, but in practice they were not much different and 
could be computed with significantly less effort; they were 
valuable both as diagnostic tools and as components of more 
complicated analyses. We found that the marginalized like- 
lihoods in the two-band EE /BB case were not significantly 
different from the conditional likelihoods, which was expected 
from the orthogonality of EE and BB due to the uniform par- 
allactic angle coverage of the scan strategy. The correlation 
between the EE and BB bands was less than 5%. 

We assigned errors to EB bandpowers by fixing EE and 
BB to their maximum likelihood values (computed assuming 
EB = 0) and then computing conditional likelihoods in EB, 
rather than marginalizing over EE and BB. The reason is that 
we found EB bandpowers to be correlated significantly with 
the others, but there was a strong prior belief that EB = 0. 
Likewise, when we assigned errors to EE and BB bandpowers, 
we did not include additional uncertainty from marginalizing 
over EB. 

For the 14-band analysis, we sampled the marginal distribu- 
tion using Markov Chain Monte Carlo (MCMC). The specifics 
differed for the two pipelines and are discussed briefly in the 
following; further details are given in Smith ( 2007| l and Hyatt 
( |2008l l. 

5.3.1. Pipeline 1 

Since the signal and noise covariance matrices were dense 
matrices of size n « (2 x 10"*), and likelihood operations typi- 
cally involve 0{n^) matrix operations, several optimizations 
were required to make the MCMC analysis computationally 
feasible. 

The first optimization was the use of signal-to-noise eigen- 
modes to reduce the effective matrix size; after a few 0{n^) op- 
erations to pre-compute the eigenmodes and signal covariance 
matrices, subsequent likelihood operations could be performed 
using matrices of smaller size <C n. For CAPMAP, we 
found that no information was lost by choosing n^^ = 4500, 
reducing the dimensionality by more than a factor of 4. 

The efficient calculation of conditional likelihoods enabled 
a fast MCMC technique using Gibbs sampling. The joint like- 
lihood £(ai, ...,Q!m) was sampled by sampling the conditional 
likelihoods repeatedly. After the nth iteration, we had a new 
MCMC sample (a*"', • • • , a^') from the joint likelihood, and 
we saved the conditional likelihoods C(ai\a^j'^j). One can show 
that the conditional likelihood averaged over MCMC samples. 



''sampler 



samples 



(14) 
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converges to the marginalized likelihood C{aj). By generating 
a full conditional likelihood in each iteration, rather than just 
generating samples, the rate of convergence was dramatically 
improved. We found that only w 10^ MCMC samples were 
needed to obtain convergence. 

5.3.2. Pipeline 2 

We used the standard Metropolis-Hastings algorithm 
( Metropolis et al.|1953[ Hastings] 1970[ l with a Gaussian pro- 
posal density. We matched the proposal density closely to the 
actual marginal distributions by first calculating the conditional 
likelihoods, obtaining acceptance probabilities as high as 70%. 
We used 350000 MCMC samples to obtain the final results. 
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although convergence was typically achieved with 3-5 times 
fewer samples. 

Because of the Fourier space optimization, it was also feasi- 
ble to calculate the joint likelihood by a brute force evaluation 
of equation ^T2\ on a 60*'-element grid for M < 4. This en- 
abled us to verify that our MCMC implementation correctly 
reconstructed the marginalized likelihoods. 

5.4. Simulations 

All aspects of both pipelines were tested extensively on 
simulated data sets. To create ?i full-season simulation, we pro- 
cessed the complete data set, replacing each 100-Hz sample 
for all 44 channels and all 280000 cycles with a data point 
drawn at the same RA, declination and parallactic angle from a 
simulated CMB realization of the CAPMAP field, with Gaus- 
sian noise added at the same level present in the real data at 
that time. The simulated timestreams were then processed into 
simulated map vectors by the binning and mapmaking codes, 
applying the same data selection used for the real data. The 
resulting map vector had the same white noise properties as 
the actual data set, but a known underlying CMB signal, which 
could be set to have any desired power spectrum, typically the 
concordance model for EE and zero power for BB. 

These simulations enabled end-to-end tests of the binning, 
noise estimation, mapmaking, and power spectrum estimation 
schemes. The maximum likelihood estimator algorithms were 
verified to be unbiased at the subpercent level (much smaller 
than the statistical error bars), and the asymmetric confidence 
intervals were verified to have the expected statistical proper- 
ties, in all I bands. Timestream simulations were exchanged 
between the two pipelines, verifying that they were consistent 
with each other as well as internally. 

For pipeline 2, it was also convenient to generate simulated 
receiver maps in the {6, LST) binning space directly for many 
studies. Each such simulation could be generated in only a few 
minutes, and the mapmaking and power spectrum estimation 
proceeded with similar speed in the half-Fourier space. 

5.5. Null Tests 

The pipelines were designed to allow a suite of null tests — 
statistics that explored a large set of possible sources of con- 
tamination and systematic uncertainty. The first step was to 
make separate maps (A^i ,Xi) and (A^2,X2) for two disjoint sub- 
sets of the data (where A^,- denote the noise covariance matrices 
and X, denote the map vectors). We then constructed a null 
map (A^,x) as: 



N = Ni +N2 

X = Xi-X2, 



(15) 



and we performed a two-band power spectrum analysis on the 
null map, fitting both EE and BB as flat bandpowers. The null 
maps were constructed to have no contribution from the CMB 
and thus be consistent with zero signal. This procedure tested 
for systematic contamination associated with the way in which 
the data were divided into subsets, e.g. sun contamination in 
the case of a day /night split. Furthermore, the null maps could 
be examined at early stages of the analysis without biasing the 
final results. 

The differencing procedure defined by equation ( fTS) is ap- 
propriate when the two maps have the same spatial resolution. 
This was the case for almost all of our null maps, but there 
was one important exception: the Q-W null map, which tests 
for overall agreement between the two frequency bands. In 



this case, we handled the distinct beam sizes by degrading 
the W-band map to the coarser Q-band beam size before dif- 
ferencing using equation ( [TS] ). The degrading operation was 
performed directly in map space and applied consistently to 
the noise covariance matrix as well as the map vector. Since 
the null map was limited already by the coarser resolution, 
no information was lost in this procedure. Map differencing 
with distinct beam sizes was only implemented in pipeline 1 , 
because this procedure mixes Fourier modes. 

We summarized likelihoods resulting from null maps with 
an equivalent , which was obtained by calculating the prob- 
ability to exceed zero and expressing it as a corresponding 
with one degree of freedom. 

6. POWER SPECTRA 

In this section we present the power spectrum estimates 
for EE and BB. Power spectrum analysis of the Q-W dif- 
ference map and fits for EB power are also presented along 
with results of the comprehensive null test suite. Our primary 
measurements are of the EE and BB power spectra in the range 
200 < ^ < 3000. 

Q+W Spectra — In Figure |9] we present results from each 
pipeline for a likelihood analysis of the Q- and W-band maps 
together Although the two pipelines were independently devel- 
oped and implemented, and differed in key details (including 
pixel size and shape, map space, and noise modeling) as de- 
scribed in ^ the two sets of results agree very well. The 
small variations between the two are consistent with expected 
statistical variations from the different binning and mapmaking 
procedures, as verified by full-season simulations exchanged 
between the two pipelines. 

The results also agree with the concordance model: the EE 
power is significantly detected in multiple bands and BB power 
is consistent with zero. The window functions for the seven 
bands have been calculated as in (Knox ( 1999) and are plotted 
in panels under the power spectra in Figure [9] The effective I 
for each band is taken as the maximum of the window function, 
and the horizontal error bars indicate the 68% highest-density 
confidence interval. At high I, almost all the statistical weight 
came from W-band because of its smaller beam size; at low i, 
the weight was slightly dominated by Q-band because of its 
larger survey area. 

Two-band Spectra — Figure [TO] shows the joint likelihood con- 
tours obtained by one of the pipelines for the EEIBB two- 
band analysis. The EE band yielded an 1 Icr detection with a 
concordance multiplier 1.36;!;q 22- obtained a flat 

bandpower 0.6^j 5 pK^, with a 95% confidence upper limit 
4.8 pK^. 

Q- and W-band Spectra — In Figure [TT| we show power spec- 
trum estimates obtained by analyzing Q- and W-band maps 
separately. The results are in good agreement with each other 
and with the combined Q+W power spectrum shown previ- 
ously in Figure [9] providing evidence for a lack of foreground 
power 

Q-W Spectra — We next consider a map obtain ed by differenc- 
ing the Q- and W-band maps, as described in §5.5 Since the 
calibrations for each constituent map included the frequency- 
dependent thermodynamic correction (see ^4.6| l appropriate 
for measuring fluctuations from a 2.7 K blackbody, this map 
is insensitive to CMB fluctuations but would be affected by 
foregrounds. Moreover, when power spectra are estimated 
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3000 



Fig. 9. — Q+W combined results from pipeline 1 (squares) and pipeline 2 (circles) for EE (left) and BB (right) fiat bandpowers. The points for pipeline 1 (2) are 
shifted by £ = —10 (10) for clarity. The EE concordance model is show n as t he solid curve, and the window functions for each band are plotted below the power 
spectrum results. Neighboring bands have only small correlations (see ^5.3^ . The error bars show the statistical errors only and do not reflect the 20% overall 
calibration uncertainty or the small systematic eiTors discussed in ^ 
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Fig. 10. — EEIBB joint likelihood contours from pipeline 2 for the 
two-band Q+W fit. Contours are drawn at the one-dimensional na levels: 
Cp exp(-;i^/2), where Cp is the peak likelihood. 

from this map, their errors do not include sample variance 
contributions from the CMB. Thus, analysis of the Q-W map 
probes for foreground contamination more deeply than does 
checking for consistency between the power spectra of the W- 
and Q-band individual maps. 

Figure 12 presents EE and BB estimates from the Q-W map. 
No significant signal is apparent. We evaluated the hypothesis 
that all the EE {BB) bandpowers were consistent with zero 
by summing their equivalent statistics; we found 7.4 (8.7) 
for the 7 degrees of freedom. The lack of power in the Q-W 
map suggests that the polarized foregrounds at 40 and 90 GHz 
are smaller than the noise level of the null map (w 90 /iK- 
arcmin). The null result also confirms that the maps from the 
two frequencies can be reliably combined to produce the Q-l-W 
spectra in Figure [9] 

EB Spectra — A nonzero EB signal would violate parity in- 
variance and is therefore expected to be absent cosmologically 
and to be generated only by a limited set of systematic effects, 
such as an overall error in the detector polarization angles; of 
course, foregrounds can also generate an EB signal. (In con- 
trast, BB is parity-invariant and can be generated not only by 
cosmological effects such as gravity waves and gravitational 
lensing, but also by a wider range of systematics.) Figure [13] 
shows estimates of seven bandpowers for the EB cross power 
spectrum. The to zero is 6.2 for seven degrees of freedom, 
verifying the expected absence of signal. 



The Null Test Suite — We evaluated a set of null tests as de- 
scribed in %5.5\ each based on dividing either the W- or the 
Q-band data into disjoint subsets to form null maps, which fell 
into the following five classes: 

1 . Time divisions: difference maps among the three peri- 
ods, between day/night data, between morning/evening 
data; and quadrature data. 

2. Cryostat divisions: data from all radiometers in one 
cryostat versus all the others. 

3. Spectral divisions: one frequency subband versus the 
others. 

4. Housekeeping divisions: differences between the halves 
of the data with the highest/lowest IF amplifier, LNA, 
and telescope mirror temperatures. 

5. Performance divisions: maps from the polarimetry chan- 
nels with the best/worst performance in: 1// noise, 
ground-synchronous structure, and I ^ Q leakage. 

In total there were 72 null tests performed on the data: 17 
for Q band, 19 for W band, each for EE and BB. Figure [T4] 
displays the bandpower results from two-band fits for all of 
the tests; they scatter around zero power To assess the perfor- 
mance of the null test suite, for each of the 72 likelihoods we 
calculated the equivalent y^. The distribution of this quantity 
is plotted in Figure [15] it closely follows the expected distri- 
bution for the hypothesis of zero power throughout the null 
test suite. The null tests were not fully independent; correla- 
tions do not bias the total but do increase its variance. The 
null test suite was studied in full-season simulations and the 
distribution of the resulting total in fact agreed well with 
the expected distribution for a x^ with 72 degrees of freedom, 
indicating that the correlations did not significantly affect the 
interpretation of this statistic. We computed the total x^ for the 
null suite as 77.8 (the corresponding probability to exceed is 
30%), in confirmation of the hypothesis that the null maps are 
uncontaminated. The summed x^ for subsets of the tests (W- 
or Q-band alone, just EE or BB) revealed no problems, nor did 
further null tests from the Q-l-W data or from EB bandpowers. 

Further Tests — Several further tests were performed to verify 
the robustness of the results. We investigated the stability of the 
data under ring mode removal by analyzing with 3 and 7 modes 
removed, rather than 5, and by turning off ground-synchronous 
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3000 



Fig. 1 1 . — W-band (top) and Q-band (bottom) power spectra from pipeline 1 (squares) and pipeline 2 (circles). EE is on the left and BB is on the right. We have 
omitted the 2000 <t< 3000 bandpower in Q-band, where the resolution is not adequate to estimate power on such small scales. 
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Fig. 12. — Power spectra from pipeline 1 of the Q— W map. The lack of power in both EE (left) and BB (right) is a good indication that, within the noise level of 
this map, polarized foregrounds at either frequency are not contributing. 
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Fig . 13. — Estimated EB bandpowers from pipeline 1 for the Q+W combined 
map, showing consistency with EB = 0. 

Structure removal. In all three cases, the results changed at 
less than the 0.1 cr level. We conclude that the data suffer 
negligible contamination from scan-synchronous structure, a 
validation of the efforts to reduce sidelobes described in ^ 
We similarly analyzed the data with a variety of different cut 
masks, corresponding to various choices of the data selection 
parameters, and found that the result varied negligibly and 
within expectations from the statistics of the samples. These 
expectations were calibrated by following the same procedure 
on a set of full-season simulations. Finally, we additionally 
analyzed the data using the alternative responsivity model 



(model 1 as described in 5 4.7 1 and found the results changed 
only by a small amount in agreement with the systematic errors 
quoted in ^ 

From the checks and null tests described so far, our results 
show no evidence of contamination from either instrumental 
systematics or CMB foregrounds. To make this quantitative, in 
the following two sections we describe our modeling of possi- 
ble instrumental systematics and of foreground contamination. 

7. SYSTEMATIC UNCERTAINTIES 

Extraction of power spectra from the CAPMAP raw data 
required the use of several imperfectly determined calibration 
parameter sets, including 1) the responsivities, 2) the beam 
sizes, 3) the polarization angles, and 4) the pointing model 
parameters. We outline in ^how we determined our best 
estimates of these parameters, which we call the baseline pa- 
rameters. Below, we describe how simulations were used to 
estimate the systematic effects on the power spectrum caused 
by measurement and modeling uncertainties for each of these 
parameter sets, and we then describe a few additional system- 
atic effects that were investigated but found to be negligible. 

The analysis codes (for mapmaking and power spectrum 
estimation) took the calibration parameters as inputs; the re- 
sults from ^ are derived using the baseline parameters. Of 
course, the simulations also required calibration parameters 
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Fig. 14. — Results from two-band fits to the suite of 34 null maps for Q-band (left) and 38 for W-band (tight). Flat bandpowers with Icr and la error bars are 
shown for each. Each null test is described further in the text. 

variations. To evaluate the impact on the power spectrum of 
neglecting the faster variations, we generated 100 full-season 
simulations and, for each, compared results from the standard 
analysis to results from analysis using the alternate respon- 
sivity model, which included short-timescale variations. The 
resulting systematic error estimate for each bandpower was 
typically 10-20% of the statistical error, increasing the overall 
error in the worst case by 3% after adding in quadrature. 

To bound the effects of these inconstant responsivities, we 
also analyzed the simulations under the unreasonable assump- 
tion that the responsivities were constant for the duration of 
the season. In this case — in which we did not account for the 
period 11 anomaly, the time- varying opacity, or even for the 
well-understood gradual responsivity decrease over the course 
of the season — the derived systematic error only doubled. 

The Q+W result was also subject to an additional systematic 
error (not affecting Q- or W-band individually) due to the 
relative calibration uncertainty between the two frequencies, 
which was dominated by the uncertainty on the W-band beam 

This was studied in receiver 
the relative calibration within the 
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Fig. 15. — The distribution for fitting each test in the CAPMAP null test 
suite. The smooth curve shows the theoretical distribution for one degree of 
freedom. The agreement indicates that the ensemble of null tests shows no 
evidence for contamination. 

to translate CMB realizations into detector timestreams or re- 
ceiver maps. Our usual approach was to generate hundreds of 
simulations using the baseline parameters, and then to analyze 
each one using perturbed estimates of the calibration param- 
eters in the parameter set under study. These new estimates 
were drawn from distributions consistent with measurement 
and/or modeling uncertainties. In addition, each simulation 
was also subject to the standard analysis. The two analyses of 
each simulation were compared, and the shifts in the likelihood 
peaks for each of the 14 £ bands were recorded. 

For every parameter set studied, the mean peak shift was 
consistent with zero. The standard deviations of the peak shift 
distributions were used to estimate the systematic uncertain- 
ties, which are listed in Table[T|for the Q+W EE bands. We 
also estimated systematic errors for W- and Q-band separately, 
finding similar results. The table does not list the BB system- 
atic effect levels; they are comparable, in comparison with the 
statistical error bars, to the EE levels. 



efficiency discussed in 5 4.6 



map simulations by varying th 
w ±5% distribution expected from the optical simulations and 
was found to be subdominant in all cases to the responsivity 
model uncertainty. The first column of Table [T] reports the two 
contributions added in quadrature. 

These systematic error estimates were perhaps smaller than 
might have been predicted in light of the 10%-level effects 
ignored in the baseline responsivity model. Two features of 
CAPMAP provided immunity: the large number of channels, 
and the fact that the entire sky was surveyed every six hours, 
which diminished the effect of short-timescale variations. We 
note that though these systematic errors are small relative to our 
statistical errors, they are not small compared to predictions of 
the primordial BB spectmm; fortunately this systematic effect 
would be greatly reduced with successful temperature control 
to stabilize the responsivities. 



7.1. Responsivity 

The time dependence of the responsivity for each polarime- 
ter was modeled as described in |4.7| We considered two mod- 
els and chose the simpler one, which ignored short-timescale 



7.2. Beams 

Observations of Jupiter were used to find the beamsize for 
every receiver to a few percent, as described in ^4.2| Each 
W-band (Q-band) beamsize was within 5-6% of the W-band 
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(Q-band) mean. We took advantage of the good agreement 
by treating all the W-band (Q-band) channels as having the 
same beamsize, which facilitated combination of the data from 
separate polarimeters in the mapmaking stage. We evaluated 
the systematic impact of this approximation by generating sim- 
ulations in which each receiver map resulted from convolution 
with the actual measured beamsize of that receiver, and then 
subjecting the maps to the standard analysis. The resulting 
systematic error was not the dominant effect in any t band, 
validating our simplified treatment of the beams in the pipeline. 

We also investigated a possible uncertainty from ignoring the 
< 8% elongations of the beams. We numerically calculated the 
window function (before the complication of scanning) for the 
worst elongation, and then differenced it from the symmetric 
one used in the analysis. Properly folding this difference with 
the concordance model EE power spectrum showed that, for 
all bands, the corresponding correction was a small fraction of 
the beam uncertainty discussed above so that it could be safely 
neglected. 

Finally, we used a similar method to evaluate the validity of 
our treatment of the beams as Gaussian, in view of the side- 
lobes described in ^4.6| We modeled the sidelobes by adding a 
second Gaussian to the main one, calculated a simplified win- 
dow function, and again found this effect negligible compared 
to the beam uncertainty estimated above. 

7.3. Pointing 



Results from the pointing solution are described in § 4. 1 The 



largest uncertainty was in determining the global pointing: the 
position of the center of the survey region (nominally the NCP). 
The associated systematic effect was evaluated from 500 re- 
ceiver map simulations, each analyzed with an estimate for the 
azimuth and elevation of the survey center drawn from Gaus- 
sian distributions centered at the NCP with widths 18" and 29" 
respectively. The resulting systematic error was not the domi- 
nant effect in any I band. The offsets of the polarimeters from 
the array center were better constrained by measurements than 
was the global pointing; simulations confirmed that they led 
to a negligible additional contribution to the power spectrum 



7.4. Polarization Angles 

The detector polarization angles were measured with a typ- 
ical error of 1.5° as described in ^4.4| We evaluated the sys- 
tematic errors reported in Table [T] with 500 receiver map simu- 
lations, analyzing each one with a set of detector angles drawn 
from Gaussian distributions with a =1.5° about their baseline 
values. We also considered the effect of a global bias in the 
detector angles, which would rotate /i-modes into B-modes 
and generate non-zero EB power The few-degree precision of 
our measurements of each detector polarization angle limited 
the possible magnitude of a global bias sufficiently to make its 
impact neghgible. 

7.5. Other Effects 

We explored a few other potential sources of systematic 
uncertainty and found them to be smaller than the others. These 
have not been included in Table [T] instead, we indicate here 
their small impacts by giving the uncertainty estimated using 
a two-band analysis of simulations for Q+W. We report the 
uncertainties as pairs: in which X is a multiplier to the 

concordance model for EE and Y is the uncertainty for BB in 



As discussed in 5 4.5 there was an additional systematic 
uncertainty from the opacity determination arising from un- 
modeled responsivity variations in the total power channels 
used in its measurement. This uncertainty was estimated by 
analyzing full-season simulations with an alternative opac- 
ity model relying on different total power channels with sig- 
nificantly different responsivity variations. The results were 
[0.005, 0.03]. 

We also studied the effects of the I ^ Q leakage described 



in 5 4.3 Worst-case estimates for the effect of the monopole 
and dipole terms predicted that they would contribute at most 
2% to measurements of the EE power spectrum (McMahon] 
2006 1. We simulated the effects of the monopole term by 
adding appropriate amounts of TT to the polarization maps 
and found the uncertainties to be [0.003, 0.04]. The EE ef- 
fect was consistent with zero, with an uncertainty that was 
nearly a factor of 10 smaller than the 2% estimate because 
of suppression from the ring scan strategy: each spot on the 
sky was measured with many parallactic angles. Quadrupole 
leakage is not suppressed by scanning or sky rotation, but 
the CAPMAP optics were specifically optimized to keep the 
quadrupole small. We estimated the quadrupole contribution 
to the measured EE spectrum to be less than 0. 1 %, smaller 
than the monopole systematic. 

Pipeline 2 was subject to an additional systematic error 
due to the assumption of LST-invariant noise ( ^5.2| l. The 
effect of ignoring the w 10% variation of the noise with LST 
was studied in receiver map simulations, which showed no 
detectable bias. The uncertainties on EE and BB were [0.004, 
0.02]. 



As mentioned in 5 5.1 some of the polarimeters suffered 



enough 1 // noise to perturb slightly the distributions of their 
statistics from the five-ring-mode fits performed on each 
cycle during the data selection process, resulting in a high 
grand x^. The analysis pipelines assumed white noise (5 5.2 1. 
We evaluated the systematic impact of this assumption by gen- 
erating full-season simulations including 1 // noise (using the 
method of Plaszczynski ( 2005| ) in each polarimeter consis- 
tent with measured noise spectra, and subjecting them to the 
standard analysis. We found a negligible bias, and the system- 
atic uncertainty was [0.003,0.05]. Finally, we note that the 
parameter-division grand null test (equivalent to differenc- 
ing maps made from channels with high and low levels of 1// 
noise) passed at the 1 a level for both EE and BB and in both 
frequency channels (Fig. 14 1. 



7.6. Systematics Summary 

The rightmost column of Table [T| gives the total systematic 
error added in quadrature, first in /iK^ and then in units of 
the statistical error. The largest systematic effect relative to 
the statistical error is for the EE 1401-1700 band. There the 
systematic error is 37% of the statistical error, increasing the 
total error for that band by 6.6%. 

We have also calculated the total systematic effect from all 
the contributions considered above (by adding the errors in 
quadrature) for the two-band fits as an overall summary: the 
result is 2% of the concordance model for EE and 0.13 /j,K^ 
for BB, both in the 2-3000 band. 

In addition to these effects, there is an overall multiplicative 
uncertainty of 20% (16%) in W-band (Q-band) on both the 
power spectrum results and their error bars, after adding in 
quadrature the 9.4% (7.6%) absolute uncertainty on the re- 
sponsivity calibration and the 2% uncertainty on the opacity 
determination (and taking the appropriate square since the re- 
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TABLE 1 
Q+W Systematic Effects 



Band 


Responsivity 


Beams 


Pointing 


Angles 


Total 


1 


0.6 


0.2 


0.3 


0.2 


0.7 (0.15(7) 


2 


1.4 


0.5 


0.7 


0.3 


1.7 (O.lScr) 


3 


1.8 


1.2 


1.3 


0.9 


2.7 (0.21(7) 


4 


2.5 


1.7 


1.9 


2.0 


4.1 (0.27(7) 


5 


4.8 


3.4 


3.9 


3.7 


8.0(0.37(7) 


6 


6.6 


4.0 


4.5 


4.8 


10.1 (0.29(7) 


7 


5.3 


3.1 


3.7 


3.0 


7.8 (0.29(7) 



Note. — Systematic error levels are given for all 7 Q+W EE power 
spectrum evaluation bands; the errors for W- and Q-band separately or for 
BB are comparable relative to the respective statistical errors. All errors 
are quoted as flat bandpowers in /xK- and are calculated as described in ^ 
The rightmost column gives the total systematic error added in quadrature; 
this quantity is to be added in quadrature with the statistical errors given 
in Table I2] Shown in parentheses is this same quantity divided by the 
half-widtn of the \a confidence interval, which allows the systematic 
uncertainty to be compared with the size of the statistical error bar. Not 
reflected in these numbers is an overall 20% calibration uncertainty 
affecting the results and the error bars equally. 

suits are expressed in units of /iK^). The overall calibration 
uncertainties on the Q+W results varied with £ and were stud- 
ied in simulations. At low £, Q- and W-band contributed equal 
weight to the result, producing a net calibration uncertainty of 
18%. At high £, only W-band contributed, and the calibration 
uncertainty is 20%. These errors do not affect the detection 
significance. 

8. FOREGROUNDS 

From results presented in ^ we find no evidence for signif- 
icant foreground contamination in our data: no power in the 
Q-W map and no significant BB power in any of our maps. 

We also tested the Q- and W-band maps for possible point 
source contamination by excising pixels in which the measured 
signal exceeded the mean by 3.5a (approximately 70 /iK at 
the 3.4' resolution used for this study). There were 4 (2) such 
pixels in the Q-band (W-band) map; masking these pixels 
produced a negligible change in the multi-band power spectra. 

Such tests have limited statistical power, however, so we 
supplemented them with simulations. Here we describe our 
methodology ( Vanderlinde||2008 1 and results for both point 
sources and diffuse galactic contributions. 

Using the full-sky map of synchrotron radiation provided by 
Haslam et al. ( 1981| l, we extrapolated with /3 = -3 to produce 
synchrotron intensity templates at 40 and 90 GHz. Similarly, 
following the method of Finkbeiner et al. ( 1999 1, which models 
thermal dust as two distinct populations with spectral indices 
of 1.7 and 2.7, dust intensity templates were produced at the 
two frequencies. Only the Q-band power from these templates 
needed to be scaled down by a factor of about four to match that 
reported by WMAP ( Kogut et al.|2007| l within the CAPMAP 
survey region. 

The polarization angle at each point on the sky was chosen 
following |Giardino et al.| ( [2002) l, and polarization levels of 10% 
and 5% were taken as conservative estimates for synchrotron 
and thermal dust, respectively. We took the two polarizations 
to be aligned, as a worst but not implausible case. 

The WMAP point source catalog ( Hinshaw et al.|[2007[ ) 
is currently the only one that includes the NCP (at galactic 
latitude 27.1°). It contains no sources within the CAPMAP 
survey region, to thresholds of roughly 300 fiK and 500 /iK in 
Q- and W-band, respectively. 

Point sources bright in microwaves and outside the galactic 
plane are dominated by Active Galactic Nuclei, the vast major- 



ity of which are blazars (Giommi & Colafrancesco, 2006): flat 
spectrum {a w -0.25, where Si, ^ v") and wea kly (< 10%) 
polarized. A recent comprehensive blazar survey ( Healey et al. 
2007| l found no sources with flux at 4.8 GHz greater than 
65 mJy (~ 200 /iK for a W-band receiver, assuming a = -0.25) 
within the CAPMAP survey region. 

To proceed, we simulated the effects of fainter sources, 
guided by the well-measured statistics of the blazar popula- 
tion (s ee e.g. [Giom mi & Colafrancesco ( 2006}, [Giommi et al.| 
( 2006^ , de Zotti et al. 1 (2004)). We assumed a polarization level 
of 3%. Our simulations agree well with |Tucci et al. (2004i, 
who estimate power spectra from extragalactic radio sources, 
taking into account measured source polarization distributions 
at 1.4 GHz, and including both steep and flat-spectrum sources. 

For each of 100 simulations, a foreground realization was 
added to a map with a CMB and noise realization matching 
the final data set. For W-band, the mean change in power 
was never more than 1 fiK^ in any band. In Q-band, the 
effect was at the level of 5% of concordance; we find, for 
example, at ^ = 1000 about 1 /iK^ from faint point sources 
and 0.5 fiK^ from synchrotron radiation. We conclude that 
our W-band results are likely foreground free while those in 
Q-band could be contaminated at tiie level of w 5%. However, 
the simulations show that our main Q+W results should be 
contaminated at less than 2%. 

9. PIPELINE COMPLEMENTARITY AND COMPARISONS WITH 
OUR PREVIOUS WORK 

In this section we summarize the roles of the two pipelines. 
We follow with summary comparisons to our previous work. 

9.1. Complementarity of Pipelines 

Our work has benefited from having two pipelines in a 
number of ways. The agreement of the results from the two 
pipelines is strong evidence that the CMB results are robust. 
The two pipelines traded simulations to verify their code im- 
plementations. As described in ^5.1| each initially developed 
independent data selection methods, though the cut variables 
developed by pipeline 1 were discovered to be the more effec- 
tive and were eventually adopted by both pipelines; it was also 
discovered that it was necessary to include the longer timescale 
explored by pipeline 2. 

Having two pipelines allowed us to subject the data and 
its analysis to a wide battery of tests; neither pipeline had to 
implement every test. 

Pipeline 1 was in a sense a conventional implementation. It 
was capable of producing likelihood estimators with arbitrary 
noise realizations and without significant approximation. How- 
ever, when the pixel count exceeded « 10'*, the processing time 
became prohibitive. The majority of studies with this pipeline 
were done with about 2500 identically sized 3.4' pixels, the 
size of the W-band beam. These studies included: the evalua- 
tion and optimization of possible scan strategies; development 
of the null-test suite, the Q-W and EB null tests, development 
and study of the selection criteria; and preliminary studies of 
the effects of 1// noise in the timestreams. 

Pipeline 2 was more novel, exploiting the speed advantages 
of the half-Fourier technique. The calculation of the likelihood 
from the data vector at one point in parameter space needed 
less than 1 CPU second. This enabled the fast simulation and 
analysis of a variety of systematic effects with 1.2' pixels, 
including those reported in Table [Tj in particular, evaluation 
of the impact of responsivity uncertainties required hundreds 
of full-season simulations from pipeline 2. Other studies with 
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Q+W Power Spectrum Results 
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Fig. 16 . — Compilation of all measurements to date of the CMB EE power spectrum. Values are taken from the pre sent work (pipelin e 2) and from MAXIPOL 
jWetal. 2007 1, DASI ( Leitch et al. 2005 1. CBI ( Sievers et al. 2007 1, BOOMERanG ( Montroy et al. 2006 1, WMAP jPage et al.|2007) . QUaD jAde et al.pOOT) , 
and the CAPMAP first season (Barkats et al. 2004 1. The error bars show statistical errors only in all cases. 

this pipeline included evaluation of the null test suite with 1 .2' 
pixels with many different data selection schemes, simulation 
of the null test suite, and study of the impact of 1 // noise 
with full-season simulations. Pipeline 2 required the noise 
for each polarimeter to be approximately independent of LST 
Although the noise condition was satisfied for the full-season 
data set and was also adequately satisfied for the three periods, 
this limitation prevented pipeline 2 from performing such tests 
as the day/night null test and the foreground study in which 
several possible point sources were excised. 

Pipeline 2 ran routinely with 1.2' pixels and determined 
14-band conditional likelihoods in 48 CPU hours; pipeline 
1, running with 1.7' pixels, took about 175 CPU hours for 
the same task (1.7' pixels were small enough, given the noise 
levels). The final results were produced with the MCMC 
method; the implementations were different, with pipeline 1 
taking about 10000 CPU hours and pipeline 2 about 1500 CPU 
hours. 

9.2. Comparisons with Our Previous Work 

Our first publication reported a 2a detection over one large 
band obtained from one season's worth of data with 4 W-band 
radiometers. This work is a significant improvement over 
that result, most importantly because we operated for more 
time with a greatly expanded array and introduced a second 
frequency. 

To capitalize on the improved statistics, for this work we 
treated systematic effects in much more detail, including in- 
corporation of the beam efficiency effect that we had taken to 
be negligible in the past; in fact the magnitude of that effect 
was equal to the calibration error we reported in Barkats et al. 



Note. — All results are in /jK^ and are expressed as maxi- 
mum likelihood flat bandpowers with 68% confidence intervals. 
The effective multipoles are listed for the EE bands; those for 
BB differ only slightly. The errors given are statistical errors 
only; the systematic errors are discussed in ^ 

systematic effects. One was the improved AR coatings on the 
lenses (^and which greatly reduced the polarimeters' 
sensitivity to unpolarized point sources (subdominant to the 
expected contribution from polarized point sources). Another 
was the reduction in telescope sidelobes (^, which made 
ground-synchronous structure (5 5.1 1 negligible, eliminating 
concerns that variations in such residual scan-synchronous sig- 
nal might contaminate the data. Finally, the highly redundant 
and symmetric ring scan strategy provided another layer of 
immunity. 



(2004 1. Several improvements reduced our susceptibility to 



10. FINAL RESULTS AND CONCLUSIONS 

Figure [16] shows our Q+W EE results together with those 
of all other experiments that have reported EE measurements. 
Both EE and BB bandpowers are listed in Table [2] The ex- 
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Fig. 17. — 95% confidence upper limits on the BB p owe r spectrum from 
CAPMAP and other recent experiments (as cited in Fig. |16| . We show limits 
(neglecting calibration uncertainties) from all published experiments that 
report limits on BB power. Limits were calculated by integrating only the 
pai1 of the likelihood coiTesponding to Cf ® > 0. For MAXIPOL and DASI, 
published upper limit values were used, while for CBI and BOOMERanG, 
upper limit s were calculated from offset lognormal approximations ( Bond 
|et al.|206o) to the lilcelihood functions, using published parameters. QUaD 
does not report upper limits. The expected signal from gravitational lensing 
(smooth curve) sliows the next target for fine-scale CMB polarization. 

cellent agreement between the pipelines means that it matters 
little which set of results we report; we have arbitrarily chosen 
those of pipeline 2. 

We detect polarization at 1 Icr, and our highest £ bandpowers, 
while not of sufficient precision for a detection at the expected 
level, nevertheless extend the community data set into these 
scales. Our results, like all those already published, are in good 
agreement with the concordance model. Figure [17] presents 
our BB results, interpreted as 95% CL exclusions, again in 
comparison to previous work. While no experiment is yet 
sufficiently close to detection of gravitational lensing, our 
sensitivity is among the best. 

We mention again our result on the Q-W map (Fig 
Unlike the suite of 72 null tests described in ^ this one did not 
have to pass; that it did pass constitutes an additional science 
result — we find that maps from both frequencies are consistent 
with the same underlying CMB realization and hence limit 
foreground contamination. 

Our Q/U coverage was very uniform, owing to the ring 
scan strategy, so that our EE and BB bandpowers were largely 
uncorrelated. The scan strategy provided other benefits, in- 
cluding continuous monitoring of the opacity and data quality. 
We exploited the symmetry of the survey region to explore the 
speed advantages of analysis in a half-Fourier space. Our two 
complementary pipelines were fully tested with end-to-end 
simulations, and a variety of systematic uncertainties were 
estimated from detailed simulations. 

We emphasize here that our approach was conservative: our 
selection criteria, in particular, were (in part) based on the 
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performance of the suite of null tests. We note that while the 
null tests might have passed or failed depending on choices 
of cuts, the final results themselves were quite insensitive to 
cutting parameters. The situation was the same with ground- 
synchronous structure ( ^5.1| l: we projected out such a mode 
in the analysis, and reported results with the associated loss in 
sensitivity, but in fact not doing so made no significant change 
to the recovered central values. Similarly, we projected out five 
modes associated with the ring scan, but analyzing with only 
three modes removed left the results statistically unchanged 
(but with smaller errors). 

It might seem surprising that the polarization results from 
CAPMAP, which took data for just a few months and observed 
from not a particularly ideal location, compare to those from 
other experiments operating with comparable numbers of de- 
tectors and observing for longer times from sites regarded as 
the best on the planet. There are a few factors that contributed 
to this: our noise was sufficiently well understood that we were 
able to use auto- (rather than only cross-) correlation power 
spectra; the stability of our residual ground-synchronous struc- 
ture allowed us to remove ground contamination with a smaller 
number of modes than experiments that use field differencing; 
and the 7 m telescope allowed us to center our sensitivity in 
multipole space where the expected signal, for both E- and 
Z?-modes, is maximal. 

The next target at fine angular scales is the detection of BB 
power from lensing; reaching this goal will be a challenge, 
requiring more detectors and even better immunity from sys- 
tematic uncertainty. 
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